CLNS 06/1985 
CLEO 06-25 

Search for invisible decays of the Y(IS') resonance 

P. Rubin, 1 C. Cawlfield, 2 B. I. Eisenstein, 2 I. Karliner, 2 D. Kim, 2 N. Lowrey, 2 P. Naik, 2 
M. Selen, 2 E. J. White, 2 J. Wiss, 2 R. E. Mitchell, 3 M. R. Shepherd, 3 D. Besson, 4 
T. K. Pedlar, 5 D. Cronin-Hennessy, 6 K. Y. Gao, 6 J. Hietala, 6 Y. Kubota, 6 T. Klein, 6 
B. W. Lang, 6 R. Poling, 6 A. W. Scott, 6 A. Smith, 6 P. Zweber, 6 S. Dobbs, 7 Z. Metreveli, 7 
K. K. Seth, 7 A. Tomaradze, 7 J. Ernst, 8 K. M. Ecklund, 9 H. Severini, 10 W. Love, 11 
V. Savinov, 11 O. Aquines, 12 Z. Li, 12 A. Lopez, 12 S. Mehrabyan, 12 H. Mendez, 12 
J. Ramirez, 12 G. S. Huang, 13 D. H. Miller, 13 V. Pavlunin, 13 B. Sanghi, 13 I. P. J. Shipsey, 13 
B. Xin, 13 G. S. Adams, 14 M. Anderson, 14 J. P. Cummings, 14 I. Danko, 14 D. Hu, 14 
O : B. Moziak, 14 J. Napolitano, 14 Q. He, 15 J. Insler, 15 H. Muramatsu, 15 C. S. Park, 15 

E. H. Thorndike, 15 F. Yang, 15 T. E. Coan, 16 Y. S. Gao, 16 M. Artuso, 17 S. Blusk, 17 
J. Butt, 17 J. Li, 17 N. Menaa, 17 R. Mountain, 17 S. Nisar, 17 K. Randrianarivony, 17 
R. Sia, 17 T. Skwarnicki, 17 S. Stone, 17 J. C. Wang, 17 K. Zhang, 17 G. Bonvicini, 18 
D. Cinabro, 18 M. Dubrovin, 18 A. Lincoln, 18 D. M. Asner, 19 K. W. Edwards, 19 
R. A. Briere, 20 T. Ferguson, 20 G. Tatishvili, 20 H. Vogel, 20 M. E. Watkins, 20 J. L. Rosner, 21 
(N| : N. E. Adam, 22 J. P. Alexander, 22 D. G. Cassel, 22 J. E. Duboscq, 22 R. Ehrlich, 22 L. Fields, 22 

R. S. Galik, 22 L. Gibbons, 22 R. Gray, 22 S. W. Gray, 22 D. L. Hartill, 22 B. K. Heltsley, 22 
l/-) . D. Hertz, 22 C. D. Jones, 22 J. Kandaswamy, 22 D. L. Kreinick, 22 V. E. Kuznetsov, 22 

H. Mahlke-Kruger, 22 P. U. E. Onyisi, 22 J. R. Patterson, 22 D. Peterson, 22 J. Pivarski, 22 
D. Riley, 22 A. Ryd, 22 A. J. Sadoff, 22 H. Schwarthoff, 22 X. Shi, 22 S. Stroiney, 22 
W. M. Sun, 22 T. Wilksen, 22 S. B. Athar, 23 R. Patel, 23 V. Potlia, 23 and J. Yelton 23 



X 



(CLEO Collaboration) 

1 George Mason University, Fairfax, Virginia 22030 



2 University of Illinois, Urbana- Champaign, Illinois 61801 
3 Indiana University, Bloomington, Indiana 47405 
^ University of Kansas, Lawrence, Kansas 66045 
^ ■ 5 Luther College, Decorah, Iowa 52101 

6 University of Minnesota, Minneapolis, Minnesota 55455 
7 Northwestern University, Evanston, Illinois 60208 
8 State University of New York at Albany, Albany, New York 12222 
9 State University of New York at Buffalo, Buffalo, NY 14260 
10 University of Oklahoma, Norman, Oklahoma 73019 

11 University of Pittsburgh, Pittsburgh, Pennsylvania 15260 

12 University of Puerto Rico, Mayaguez, Puerto Rico 00681 
13 Purdue University, West Lafayette, Indiana 47907 

1Ji Rensselaer Polytechnic Institute, Troy, New York 12180 
15 University of Rochester, Rochester, New York 14627 
16 Southern Methodist University, Dallas, Texas 75275 
17 Syracuse University, Syracuse, New York 13244 
18 Wayne State University, Detroit, Michigan 48202 
19 Carleton University, Ottawa, Ontario, Canada K1S 5B6 
20 Carnegie Mellon University, Pittsburgh, Pennsylvania 15213 
21 Enrico Fermi Institute, University of Chicago, Chicago, Illinois 60637 
22 Cornell University, Ithaca, New York 14853 



1 



23 University of Florida, Gainesville, Florida 32611 
(Dated: December 21, 2006) 

Abstract 

We present a measurement of the branching fraction of invisible T(IS') decays, using 1.2 fb _1 
of data collected at the T(2S) resonance with the CLEO III detector at CESR. After subtracting 
expected backgrounds from events that pass selection criteria for invisible T(15) decay in T(2S) — > 
vr+vr-T(15), we deduce a 90% C.L. upper limit of B[Y(1S) -> invisible] < 0.39%. 



2 



Invisible decays of quarkonia to final state particles that cannot be detected in general 
purpose particle detectors provide a window on physics beyond the Standard Model. This 
is because the only allowed invisible decay of quarkonium states in the Standard Model is 
the decay to vv via annihilation into a virtual Z° boson. The predicted branching fraction 
for T(1S) is 

B[T(15) -> vv] = 4.14 x 10" 4 x B[T(1S) -> e+e"] (1) 

or of order 10~ 5 . Although this decay mode itself is sensitive to new physics such as the 
presence of an extra Z' gauge boson or R-parity violating effects in supersymmetric theories, 
the branching fraction remains well below current experimental sensitivity 

On the other hand, most dark matter scenarios require some coupling between the Stan- 
dard Model sector and dark matter, which could significantly increase the invisible decay 
rate of heavy quarkonia 0, 0|. Although the leading dark matter candidate in different 
theoretical models (e.g., in the Minimal Supersymmetric Standard Model) is the lightest 
supersymmetric particle with a mass typically M x > 6 GeV/c 2 , there is some evidence that 
the dark matter constituents of the universe may be as light as 100 MeV/c 2 and several 
models can be constructed to accommodate such particles. Based on a model-independent 
calculation, using only the result from the Wilkinson Microwave Anisotropy Probe (WMAP) 
[ij on the relic density of the universe, the author of 0] has predicted the decay branching 
fraction of the T(1S) to a pair of dark matter particles to be 

B[T(1S) - XX] ~ 0.41% . (2) 

Kinematic factors arising from the mass of the dark matter particles or the mediator can 
either enhance or suppress this branching fraction. Invisible T decays can produce a pair of 
dark matter particles with mass less than My/2 assuming the decay is mediated by a vector 
boson. 

Upper limits have been set on the invisible decays of the 77 and rf from J ftp — > (j>rf'> decay 
by the BES Collaboration [5]. In addition, a 90% confidence level upper limit of 0.25% 
on the invisible branching fraction of T(IS') from T(35) — ► 7t + 7t~T(IS) decay has been 
reported quite recently by the Belle Collaboration [f|, which is an order of magnitude better 
than earlier limits set by the CLEO and ARGUS Collaborations [7|. 

In order to observe the invisible decays of the T(IS') resonance, we use the T(2S) — > 
tt + tt~T(1S) transition and infer the presence of the T(IS') resonance from the missing mass 
in the event. In this analysis, we use 1.2 fb _1 data collected at the peak of the T(2S) 
resonance with the CLEO III detector operating at the Cornell Electron Storage Ring, a 
symmetric e + e~ collider. The data represent about 8.7 million T(2S) resonance decays, 
18.8% of them decaying to n + n-T(lS) @- 

The CLEO III detector has excellent charged particle tracking, electromagnetic calorime- 
try, and particle identification capabilities, with solid angle coverage for charged and neutral 
particles in the polar angle range | cos#| < 0.93 (the polar angle is measured with respect 
to the beam axis). The tracking system consists of a four-layer double-sided silicon vertex 
detector and a 47-layer drift chamber j9| residing in a superconducting solenoid that pro- 
duces a 1.5 T axial field. The Csl crystal calorimeter and the muon detection system are 
the same used in the CLEO II detector [lOj. Measurement of specific ionization (dE/dx) in 
the drift chamber provides particle identification. A Ring Imaging Cherenkov detector [ll| 
outside the drift chamber complements the dE/dx measurements, but is not effective at the 
particle momenta of this analysis. 
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The CLEO III trigger 12j makes decisions based on AXIAL trigger tracks (observed in 
the inner 16 axial drift chamber layers), STEREO trigger tracks (which extend into the 
outer 31 stereo drift chamber layers), and cluster information from the calorimeter. 

When the T(1S) resonance decay is invisible, only the two soft pions from the T(2S) 
transition can trigger the data acquisition system. The only CLEO III trigger requirement 
that is sensitive to this signal is the two-track trigger, which only demands the presence of 
at least two AXIAL trigger tracks. However, the two-track trigger has been prescaled by a 
factor of 20 in this data set, i.e., events that satisfy the two-track trigger requirement are 
counted by a prescale counter and the two-track trigger bit is only set in every 20th such 
event. This prescaling was implemented during all data taking at the Upsilon resonances in 
order to reduce the trigger rate on beam-related backgrounds. 

In order to study our selection efficiency for signal and possible background contributions, 
we use Monte Carlo (MC) simulations. The signal MC sample is created with the EvtGen 



13J generator using the decay sequence T(2S) — ► it + n T(15) followed by T(IS') — > v e v e . 



The di-pion transition is generated assuming relative S'-wave between the system and the 
T(15) as well as between the two pions. Final state radiation is simulated with PHOTOS 
|14| . Decays of T(2S) and continuum (e + e~ — > qq,q = u, d, s, c) events are generated with 
a customized version of the JETSET [15] program, while non-resonant tau-pairs (e + e _ — > 
t + t~) are generated using the Koralb/Tauola |16l.ll7j] event generator. The generated events 
pass through a GEANT-based [l8[ full detector simulation and the same reconstruction 
procedure used for the data. 

The analysis strategy is straightforward. We first select an inclusive sample of T(2S) — > 
7r + 7r~Y(lS') events by finding the pion candidates. This sample represents all possible 
T(1S) decays, including invisible decays. Then, we select the exclusive subsample of these 
for which the T(1S) decays to undetected particles. Events with an T(15) can be identified 
by the recoil mass against the di-pion (equivalent to the missing mass if only the two pions 
are considered in the event), which produces a peak at the T(IS') mass above a smooth 
combinatoric background. The total number of 7r + 7r~T(lS) events detected in the inclusive 
sample is extracted using a fit to the nir recoil mass spectrum. A similar exercise is performed 
upon the subset of the inclusive sample with no visible particles aside from the two cascade 
pions. The branching fraction of invisible decays is then determined as the ratio of these 
two yields after background subtraction. 

In order to select the inclusive sample of 7r + 7r~Y(lS) events, we require the two pion 
candidates to be oppositely charged tracks with good quality fit to a helix and with polar 
angle | cos (9| < 0.93. The tracks must originate from the interaction point within 5 mm 
in the plane perpendicular to the beam axis and 5 cm along the beam axis (luminous 
volume). The two-track trigger bit must be set in the event. In order to ensure that each 
of the two pions produce an AXIAL trigger track with high probability and together satisfy 
the two-track trigger requirement, we demand that the transverse momentum of each pion 
candidate is more than 150 MeV/c. We also demand that the magnitude of the momentum 
is less than 600 MeV/c. The recoil mass against the di-pion candidate is calculated using 
the 4-momentum of the incident beams in the lab frame (p% m = + p^+), and constructing 
the 4-product 

cm PlTTT 

Tip 

cm PlTTT ) fl 1 (3) 

M rec is required to be within 30 MeV/c 2 of the nominal T(1S) mass (9.460 GeV/c 2 ). 

A kinematic fit constraining the pion candidates to have a common vertex is performed. 
We require this vertex to have a fit of reasonable quality and be in the luminous volume 
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previously described. The improved momenta of the pions from the vertex fit are used to 
re-calculate M rec , slightly improving the resolution in this variable. 

The dominant combinatoric background in the exclusive sample with invisible T(1S) 
decay arises from two-photon fusion events (e + e~ — > e + e~7*7*) when the photons radiated 
by the incoming electron/positron beams produce a lepton pair or charged pion pair: 7*7* — > 
e + e~ I fi + fi~/7T + 7T~. Because the photon emission probability peaks in the forward direction, 
the electron/positron beam particles tend to scatter at a very small angle and leave the 
interaction region undetected. Thus only the lepton or pion pair produced by the two- 
photon fusion is detected (untagged two-photon events). The two-photon system is typically 
boosted along the beam axis and therefore the total momentum of the particles produced 
by the photon fusion tends to be parallel with the beam axis (i.e., the total transverse 
momentum is close to zero). 

In order to reduce the background due to untagged two-photon fusion events, we apply 
a cut on the polar angle of the di-pion momentum (| cos9 n7T \ < 0.9). In addition, we apply 
particle identification on the pion candidates using the specific- ionization (dE/dx) in the 
drift chamber by requiring a three standard deviation consistency with the pion hypothesis. 
This latter requirement is especially effective at suppressing the remaining two-photon fusion 
events which produce electron pairs. 

We accept multiple di-pion candidates per event if they pass the di-pion selection cri- 
teria. Random combinations of two oppositely charged tracks other than the real di-pion 
from the T(2S) — > n + ir~T(lS) transition can only contribute to the smooth, slowly vary- 
ing background underneath the peak; they are subtracted when we extract the number of 
7r + 7r~Y(lS') events from the fit to the recoil mass distribution. 

The exclusive sample of invisible T(1S) decays is selected from the initial inclusive 
7r + 7r~Y(lS) sample by the requirement that there is no additional good track coming from 
the vicinity of the interaction region and no extra good shower with more than 250 MeV 
energy in the event. Tracks with similar quality criteria to the pion candidates but with 
momentum 50 < p < 5500 MeV/c are counted as "good" tracks. Showers are considered 
"good" if they are not associated with the pion candidates. Figure [1] displays, for both data 
and MC simulations, the distribution of the maximum extra shower energy and the number 
of extra tracks for events that passed the di-pion selection and the requirement on the other 
variable. 

Based on the signal MC simulation, the efficiency of the exclusive selections on invisible 
T(IS') decays that passed the inclusive di-pion selections without trigger requirement is 
98.8%. This efficiency is in good agreement with that obtained from fully reconstructed 
T(2S) — > 7t + 7t~T(1S), T(1S) — > e + e~ /fi + fi~ data and MC simulation, as demonstrated in 
Fig. [Tb-d. In such fully reconstructed events, only tracks and showers not associated with 
any of the four final state particles are counted as extra. 

Figure [2t shows the distribution of masses recoiling against the di-pion in the inclusive 
sample. This M rec spectrum is fitted with a double Gaussian signal shape plus a linear 
background with all parameters allowed to float. The resolution of the core Gaussian is 
1.04 ± 0.04 MeV/c 2 and its fractional area is (64 ±6)% of the total. The width of the wider 
Gaussian, which arises from multiple scattering of charged particles, is twice as much as that 
of the core Gaussian. The fit gives a yield of N[7t + 7t^T(1S)} = 18, 905 ± 213 ± 400 events, 
which represents the number of inclusive 7r + 7r~Y(lS) events. The systematic uncertainty 
on the number of T(15) events in the inclusive sample is determined from variations of the 
fit to the recoil mass spectrum. We tried different bin sizes, background parametrization, 
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FIG. 1: The energy of the largest extra shower (left) and the number of extra tracks (right) in 
data and simulation. An event can enter in one of these plots only if it passes the di-pion selection 
and the cut on the other variable. Plots (a) and (b) compare the inclusive data (points) to signal 
MC simulation (solid histogram), normalized to the first data bin. Plots (c) and (d) compare fully 
reconstructed T(2S) — > tt + tt~T(1S), T(1S) — > e + e _ /// + /i~ data and the MC simulations for such 
di-lepton final states (dashed histogram) and the signal process (solid histogram), both normalized 
to the data area. 



and signal shape. 

The recoil mass distribution in the exclusive sample with no extra tracks and showers 
larger than 250 MeV is displayed in Fig. [2}3. The spectrum is fitted with the same signal 
shape by fixing the parameters of the double Gaussian to be those of the inclusive fit. The 
yield obtained for invisible T(1S) decays is 116 ± 24 ±9 events. The systematic uncertainty 
is due to fit variations. 

From the study of the generic T(2S), continuum, and tau-pair MC samples we conclude 
that the smooth combinatoric background underneath the signal peak in data with only 
di-pion selection (Fig. [2^) is mainly composed of T(2S) resonance decays (75-85%) and 
continuum e + e~ — > qq (q = u,d,s,c) events (10-15%). Because the background from non- 
resonant tau-pair events (e + e _ — > t + t~) is negligible, the remaining 5-10% of the non- 
peaking background is attributed to two-photon fusion. On the other hand, the exclusive 
selection on the extra tracks and showers suppresses the combinatoric background from 
T(2S) resonance decays, continuum, and tau-pair events to a negligible level, and therefore 
the remaining combinatoric background (Fig. [2b) is almost exclusively produced by untagged 
two-photon fusion. This is confirmed by the distribution of cos 9 nn (without selection on this 
variable), which peaks along the beam axis as expected for two-photon fusion events. 

The MC simulation of T(2S) decays also demonstrates that there is an irreducible peaking 
background in the exclusive sample due to T(1S) decays into visible particles when those 
particles lie outside the detector acceptance. According to the MC simulation, about 93% of 
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FIG. 2: Distribution of the di-pion recoil mass, M rec , in the inclusive data sample passing the 
di-pion selection (a), and in the exclusive data sample, requiring no extra tracks and showers in 
the event (b). The fit function, a double Gaussian plus a linear background, is superimposed on 
the data as a solid curve. On the lower plot, the signal shape together with the expected peaking 
background from T(LS') — > £ + £~~ (£ = e,fi,r) decays (shaded histogram) are also displayed. 



this background is caused by two body decays to e + e~ and and 5% due to decays to 

t + t~ . Other T(1S) decays that could fake an invisible signal are estimated to be negligible 
in comparison to these di-lepton final states. 

In order to get a better estimate of the dominant peaking background due to T(2S) — > 
n + 7v~T(lS), T(IS') — > £ + £~ decays, we use a MC sample generated with EvtGen using 
the correct angular correlations in decays and simulating final state radiation (the latter 
of which is particularly important for leptons). The agreement between fully reconstructed 
it + -n~ H + H~ data and the corresponding MC simulation is demonstrated by Fig. [31 

Based on the MC simulation, the acceptance (e«) of the exclusive requirement on the 
extra tracks and showers not associated with the pions for -k + -k^£ + £~ events that passed 
the inclusive di-pion selection is (8.9 ± 0.3)% with electrons, (8.5 ± 0.3)% with muons, 
and (0.8 ± 0.1)% with tau leptons in the final state. We estimate the number of invisible 
T(IS') — > £ + £~ events for each lepton flavor {£ = e, /z, r) in the data using the formula 



N[n + n-T(lS)] x x e u 



(4) 
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FIG. 3: Di-pion recoil mass, di-pion invariant mass, /i + invariant mass, and the cosine of the 
polar angle of the // + candidate in fully reconstructed T(2S) — > ir + ir~T(lS), T(15) — > /x + 'fjT 
data with two-track trigger requirement (dots) and the corresponding MC simulation (shaded 
histogram) . 



where iV[7r + 7r~Y(l 1 S')] = 18,905 is the observed number of inclusive 7r + 7r~T( IS) events in 
the data, and B m = 0.0248 is the decay branching fraction T(IS') — > /U + pT [8y, and we 
assume lepton universality (B^ = B ee = B TT ). The expected backgrounds and statistical 
uncertainties from T(15*) — > e + e~,fi + fi~, and r + r _ are 41.7 ± 1.5, 39.9 ± 1.5, and 3.8 ± 0.4, 
respectively. The total invisible background due to T(IS') decaying to lepton pairs is 85±2±4 
events, where the systematic error includes uncertainties in the B^ branching fraction, MC 
efficiency, detector simulation and reconstruction code, and MC statistics. 

After subtracting the expected background due to T(IS') — > decays from the invis- 
ible signal extracted from the fit to the exclusive data sample, we obtain an invisible signal 
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yield of 

N[T{1S) -> invisible] = 31 ± 24 ± 10. (5) 

The largest systematic uncertainty in B[T(1S) — > invisible] comes from the 32% systematic 
uncertainty in Eqn. [5j There are additional systematic uncertainties associated with the 
trigger, the di-pion finding efficiency, the fit to the inclusive data sample, and the additional 
selection criteria of the exclusive sample. The combined relative systematic uncertainty in 
B is 34%. 

When we calculate the branching fraction we assume that the trigger efficiency as well 
as the di-pion finding efficiency is the same in the exclusive and inclusive samples, and thus 
cancel in the ratio. 

There are indications that more tracks and/or showers in the final state can result in the 
two-track trigger prescale counter incrementing more than once in an event. In particular, 
we measured the two-track trigger prescale factor to be 18.3 instead of the expected 20 
in inclusive 7r + 7r~Y(lS) events by taking the ratio of signal yield from the fit to the di- 
pion recoil spectrum with and without two-track trigger requirement. We cannot directly 
measure the prescale factor for our signal events with only two soft pions since there is 
no other trigger requirement that could be always satisfied by these events. Therefore, we 
assign a 10% systematic error to the branching fraction due to uncertainty in the ratio of 
the prescale between the inclusive and exclusive 7t + tt~T(1S) events. 

We have also measured the AXIAL tracking trigger efficiency for pions with > 150 
MeV/c in fully reconstructed T(2S) — > n + iv~T(lS), T(1S) — * e + e~/[i + fi~ data where the 
extra leptons can also satisfy other independent trigger requirements and we do not rely on 
the two-track trigger for detection. We found an efficiency exceeding 99.9%, and therefore, 
we assume that the systematic uncertainty is negligible due to AXIAL trigger inefficiency. 

Based on signal MC simulation, the overall efficiency of the di-pion selection criteria 
(without the two-track trigger requirement) for invisible T(1S) decays is ~ 27%. From 
comparison to MC simulation of T(25) decays, we assign a 5% systematic uncertainty due 
to the relative difference in the di-pion finding and reconstruction efficiency between exclusive 
T(IS') decays and inclusive T(1S) decays which have significantly larger track and neutral 
multiplicity. 

The systematic error associated with the selection efficiency of the exclusive cuts on the 
number of extra tracks and showers with E > 250 MeV for signal events is estimated by 
comparing the efficiency for signal MC simulation to that for fully reconstructed T(2S) — > 
n + iv~T(lS), T(IS') — > e + e~/fi + fi~ data and the corresponding MC simulation (Fig. [Tb-d). 
Based on the differences we assign a conservative systematic error of 1%. 

After correcting the number of invisible signal events with the efficiency of the exclusive 
selection on the extra tracks and showers in the event (0.988), and normalizing with the 
number of total T(IS') decays observed in our sample of 7r + 7r _ T(lS') events (18,905), we 
obtain the branching fraction of invisible Y(IS') decays as 

B[T(1S) -> invisible] = (0.16 ± 0.13 ± 0.05)% . (6) 

The first error is statistical, from the 77% relative statistical error in Eqn. [5l and the second 
is systematic. 

We then calculate a 90% confidence level upper limit from the combined (in quadrature) 
statistical and systematic error using Feldman and Cousins' approach [l9( assuming Gaussian 
distributed error and excluding the non-physical region below zero: 

B[?(1S) -> invisible] < 0.39% (90% C.L.). (7) 
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Our upper limit is about 50% larger than the upper limit reported by the Belle experiment 
0. This is consistent with our lower statistics, which is significantly reduced by the trigger 
prescale on the two-track trigger requirement. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. D. Cronin-Hennessy and A. Ryd thank the A. P. Sloan 
Foundation. This work was supported by the National Science Foundation, the U.S. De- 
partment of Energy, and the Natural Sciences and Engineering Research Council of Canada. 



[1] L. N. Chang, O. Lebedev, and J. N. Ng, Phys. Lett. B 441, 419 (1998). 
[2] B. McElrath, Phys. Rev. D 72, 103508 (2005). 
[3] P Fayet, Phys. Rev. D 74, 054034 (2006). 

[4] D. N. Spergel et al. (WMAP Collaboration), Astrophys. J. Suppl. Ser. 148, 175 (2003). 
[5] M. Ablikim et al. (BES Collaboration), Phys. Rev. Lett. 97, 202002 (2006). 
[6] O. Tajima et al. (Belle Collaboration), |hep-ex/0611041[ 

[7] D. Besson et al. (CLEO Collaboration), Phys. Rev. D30, 1433 (1984); H. Albrecht et al 

(ARGUS Collaboration), Phys. Lett. B 179, 403 (1986). 
[8] W.-M. Yao et al. (Particle Data Group), J. Phys. G: Nucl. Part. Phys. 33, 1 (2006). 
[9] D. Peterson et al, Nucl. Instrum. Methods Phys. Res., Sect. A 478, 142 (2002). 
[10] Y. Kubota et al. (CLEO Collaboration), Nucl. Instrum. Methods Phys. Res., Sect. A 320, 66 

(1992); D. Bortoletto et al, Nucl. Instrum. Methods Phys. Res., Sect. A 320, 114 (1992). 
[11] M. Artuso et al, Nucl. Instrum. Methods Phys. Res., Sect A 554, 147 (2005). 
[12] G. D. Gollin, J. A. Ernst, J. W. Williams, and M. J. Haney, IEEE Trans. Nucl. Sci. 48, 547 
(2001); R. M. Hans, C. L. Plager, M. A. Selen and M. J. Haney, IEEE Trans. Nucl. Sci. 48, 
552 (2001); M. A. Selen, R. M. Hans and M. J. Haney, IEEE Trans. Nucl. Sci. 48, 562 (2001). 
[13] D. J. Lange, Nucl. Instrum. Methods Phys. Res., Sect. A 462, 152 (2001). 
[14] E. Barberio, B. van Eijk, Z. Was, Comput. Phys. Commun. 66, 115 (1991); E. Barberio and 

Z. Was, Comput. Phys. Commun. 79, 291 (1994). 
[15] T. Sjostrand, Comput. Phys. Commun. 82, 74 (1994). 

[16] S. Jadach, Z. Was, Acta Phys. Polon. B 15, 1151 (1984); ibid. B 16, 483(E) (1985); S. Jadach, 
Z. Was, Comput. Phys. Commun. 36, 191 (1985); ibid. 64, 267 (1991); ibid. 85, 453 (1995). 

[17] S. Jadach, J. H. Kuhn, and Z. Was, Comput. Phys. Commun. 64, 275 (1991); M. Jezabek, 
Z. Was, S. Jadach, J. H. Kuhn, ibid. 70, 69, (1992); S. Jadach, Z. Was, ibid. 76, 361 (1993); 
R. Decker, E. Mirkes, R. Sauer, Z. Was, Z. Phys. C58, 445 (1993). 

[18] R. Brun et al, GEANT 3.21, CERN Program Library Long Writeup W5013 (1993), unpub- 
lished. 

[19] G. J. Feldman and R. D. Cousins, Phys. Rev. D57, 3873 (1998). 



10 



